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Abstract 



We evaluate the Wald Noether charge entropy for a black hole in generalized theories 

> 

■^j- ■ of gravity. Expanding the Lagrangian to second order in gravitational perturbations, 

Q\ '. 

we show that contributions to the entropy density originate only from the coefficients 
and to show that arbitrary powers of matter fields and their symmetrized covariant 



of two-derivative terms. The same considerations are extended to include matter fields 



derivatives cannot contribute to the entropy density. We also explain how to use the 
linearized gravitational field equation rather than quadratic actions to obtain the same 



results. Several explicit examples are presented that allow us to clarify subtle points in 
the derivation and application of our method. 



1 



1 Introduction 



The Bekenstein-Hawking law [TJ [2] , relates the entropy Sbh to the horizon area A in units 
of Newton's constant for a black hole in Einstein's theory of gravity, 

Sbh = -rpr- • (1) 

This relation suggests that the entropy S should be purely geometric, defined strictly at the 
black hole horizon and should satisfy the first law of black hole mechanics, 

T H dS = dM . (2) 

Here, M is the conserved or ADM mass of the black hole (other conserved charges have been 
neglected for simplicity) and Tjj = k/2it is the Hawking temperature in terms of the surface 
gravity k. 

The mass M and k are well defined for a stationary black hole in any theory of gravity, 
and so their definitions do not need to be modified. Wald [31 H] proposed a definition of the 
entropy that fulfills all of the above requirements for general theories of gravity. The Wald 
entropy Sw has a clear geometric interpretation through its identification with the Noether 
charge for spacetime diffeomorphisms. Further, Sw can always be cast as a closed integral 
over a cross-section of the horizon %, 

Sw = f s w dA , (3) 
Jh 

with sw being the entropy per unit of horizon cross-sectional area. For a D-dimensional 
spacetime with metric ds 2 = gudt 2 + g rr dr 2 + J2fj~=i &ijdxidxj , dA = y/adxi . . . dx^-i ■ 

The actual Wald formula and how it comes about is briefly reviewed in Section 2. Mean- 
while, it has since been shown by three of the current authors [S] how Wald's entropy density 
Sw could also be extracted directly from the gravitational action. This, through a process 
of expanding the Lagrangian to quadratic order in the perturbations of the metric and eval- 
uating it at the horizon. The density sw can then be identified with the coefficient of the 



kinetic terms for the r, t-polarized gravitons h rt . This coefficient measures the strength of 
the gravitational coupling for the same gravitons. The advantage of this method is that it 
identifies in a straightforward way the correct units in which the area of the horizon should 
be measured to give the correct value of the entropy. It also provides a way to decide which 
terms in the expanded action can contribute to the entropy. 

Hence, one should always be able to obtain the Wald entropy by, first, expanding the 
Lagrangian around the background solution and, then, reading off the horizon value of the 
coefficient of the relevant kinetic terms; for instance, V a h rt V a h rt . Equivalently and generally 
easier to implement, one may read the same coefficients off of terms like V a V a h r t in the 
linearized field equation. 

For two- derivative theories of gravity, for instance, Einstein's and F(1Z) theories, the 
procedure is as straightforward as just described. On the other hand, for theories with four 
or more derivatives, this seemingly simple process can become rather subtle. Our current 
motivations are to provide a well-defined prescription for identifying the kinetic terms and to 
better understand why Wald's formulation still works for these higher-derivative cases. 

In this paper, we establish by explicit calculations that, for a completely generic theory 
of gravity, the kinetic contributions are indeed the only contributions to the Wald entropy. 
Further, we verify that the coefficients of these kinetic terms are always sufficient to reproduce 
Wald's formula for any number of derivatives that may appear explicitly or implicitly in 
Lagrangian. In the process we clarify the detailed properties of the Wald entropy that lead 
to such results. We then show that adding matter fields does not alter any of our results. We 
identify the correct form of the gravitational field equation that is suitable for calculations of 
the entropy via the linearized field equations. 

The rest of the paper proceeds as follows. The next section briefly reviews the standard 
derivation of the Wald formula. The new material begins in Section 3, where we consider a 
general theory of gravity and expand the Lagrangian to second order in perturbations. Then, 
using basic properties of the metric, Riemann tensor and horizon generators, we distinguish 
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between terms that can possibly contribute to the entropy and those that cannot. We go on 
to verify that the surviving terms do indeed lead to an entropy in agreement with Wald's 
expression. In Section 4, the previous considerations are then extended to generic theories of 
both gravity and matter. Our attention turns, in Section 5, to the gravitational field equation 
for a generalized theory. In Section 6, some of specific models are used to illustrate our earlier 
analysis. Section 7 contains a summary and some concluding comments. 

2 A brief review of the Wald Noether charge entropy 

In this section, we recall the derivation of the Wald entropy, following Jacobson, Kang and 
Myers j6]. Our goal is to make the paper self-contained. For brevity, we skip over the many 
subtleties and caveats in the derivations. 

One starts by varying a given Lagrangian density L with respect to all the fields {ip}, 
including the metric. In condensed notation (with all tensor indices suppressed), 



where 8 = are the equations of motion and the dot represents a summation over all fields 
and contractions of tensor indices. Also, d denotes a total derivative, so that 9 is a boundary 
term. 

Let £% be a Lie derivative acting along some vector field £. Then, given the diffeomorphism 
invariance of the theory, 5^ip = £^ip and 5^L = £^L = d(£ ■ L) . These and Eq. (0| can be 
used to identify the associated Noether current J^, 



The point being that dJ^ = when 8 = 0, and so there must be an associated "potential" 
such that = dQ^ . Now, if D is the dimension of the spacetime and S is a D — 1 
hypersurface with a. D — 2 spacelike boundary dS, then 



5L = £ ■ 5ip + d6 (<5V) 



(4) 



J 5 = e(£&)-£.L. 



(5) 




(6) 
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is the associated Noether charge. 

Wald showed [3] and later proved rigorously [I] that the black hole first law (121) is satisfied 
when the entropy is defined in terms of a specific Noether charge. Choosing the surface S as 
the horizon T-L and the vector field £ as the horizon Killing vector x (with its surface gravity 
normalized to unity [j), Wald identified the entropy as 

S w = 2tt I Q x . (7) 

*/ "hi 

Since x = on the horizon, the right-most term in Eq. (J5j) does not contribute to Sw 

One of the main advantages of Wald's formula is the simplicity of Q. To understand this, 
let us start with the Killing identity [7] 

v c VaXb = -n abcdX d , (8) 

where 1Z a bcd is the Riemann tensor. After repeated applications of this relation, the most 
general form of the integrand in Eq. (J7J) can be expressed as 

Q ab e ab = [B ab cX c + C ab cd V c X d ] e ab , (9) 

where B a b c and C a bcd are theory-dependent background tensors, while t ab = V a Xb is the 
binormal vector for the horizon. For future reference, 

tab = -tba (10) 

From the definition of e a b and since x a = (on H), it follows that the integrand simplifies 

to 

Q ab e ab = C abcd e ab e cd (11) 

or 

S w = 2n<f C abcd t ab t cd dA . (12) 



lr This particular normalization for the Killing vector will be assumed throughout the paper. 

2 More accurately, at the bifurcation surface of the horizon. This is one of the many caveats that are dealt 

with in [6]. 



We continue to sketch the analysis of [6] for theories with a Lagrangian C = £[g a b, Habcd] ■ 
First, the variation of the density L = yf—gC is found to yield 

SL = -2V a (X abcd V c 5g bd V=g) + ■■■ , (13) 
where dots are meant as terms that end up being irrelevant to the Wald entropy and 

X ahcd ee —- . (14) 

Ol-Cabcd 

It follows from Eq. (0} that Eq. (JTSJ leads to 

6 = -2n a X abcd V c 8g bd ^f + • ■ • , (15) 

where n a is the unit normal vector and 7 a & is the induced metric for the chosen surface S. 

For an arbitrary diffeomorphism 5g a b = V £j, + Vf,^ a , the associated Noether current is 
then equal to 

J = -2V a (x abcd V c {Vbia + V a tb)n a Vh) + ••• . (16) 

Let us now specialize to the horizon S — » H and (normalized) Killing vector ^ a — > \ a \ 
so that n a \/h — > e a yfa with e a = e a bX b ■ Then, using the symmetries of X ahcd (inherited 
from H abcd ^ along with Eq. (TlO|) . one can eventually translate Eq. ( fl6|) into 

J = -2W b {X ahcd V cXd t a ^) + ••• . (17) 

In this form, the potential is 

Q = _x abcd e ah e cd ^ + • • • , (18) 

and so 

S w = -2vr I X abcd e ab e cd dA , (19) 

with X abcd defined in Eq. ([HD. 
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3 Evaluating the Wald entropy 

The main goal of this section is to establish our claims that only the kinetic terms for the h rt 
gravitons can contribute to the Wald formula and that the entropy can be deduced simply 
by reading off their coefficients. We accomplish this for a generic theory by, first, identifying 
all possible kinetic terms in the quadratic expansion of the Lagrangian and, then, verifying 
that the coefficients of the relevant terms produce a result that agrees with Wald's expression. 
Along the way, we find that the contributions to the Wald entropy arise from a specific class 
of terms in the quadratically expanded Lagrangian, those terms that contain a second-order 
expansion of the Riemann tensor. 

3.1 Preliminaries 

We will begin with a pure gravitational theory. The Lagrangian for such a theory can be ex- 
pressed in terms of the metric, the Riemann tensor and its symmetrized covariant derivatives, 

£ C' \9 abi^R- abed i ^ a\R- abed (a\^ a2 

)Tl ahcd ,..] , (20) 

where the ellipsis denote increasing numbers of symmetrized covariant derivatives acting on 
the Riemann tensor. The derivatives can be expressed in such a symmetrized form, as any 
anti-symmetric combination can be converted into a Riemann tensor. 

Our objective is to expand the Lagrangian density y/^gC to second order in the metric 
perturbations, h a b = g a b ~ 9 a °b > an d then isolate the two-derivative terms which we call 
"kinetic terms" . We should also consider terms that have four or more derivatives, as clarified 
below. We will show that any kinetic term on the horizon can be expressed as 

[A abcd ] (0) Th ab V e h cd . (21) 

Here and in what follows, A aiCl2 '" represents an arbitrary tensor built out of the Riemann 
tensor, its symmetrized derivatives and the metric. A numeric superscript on a tensor denotes 
its order in fa's and V is a zeroth-order covariant derivative. 
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We can prove that a kinetic term is always of the form (|2Til . which is a four-index tensor 
with the two derivatives contracted with each other. Let us begin with the most general 
term carrying exactly two gravitons and exactly two derivatives. Repeatedly integrating by 
parts and discarding surface terms and "mass terms" which have no derivatives acting on a 
graviton, we eventually arrive at 



j^abcdef 



(0). 



(22) 



Here, the background tensor contains no explicit derivatives and is denoted by A. 



Given our Lagrangian, the background tensor 



J^abcdef 



(0) 



is built out of the tensors [g ab ] 



obi (°) 



and 



[R abcd ] {0) oc [g ac g bd - g ad g bc ] 



adbcl (°) 



(23) 



The last expression follows from the fact that, when evaluated on a stationary horizon, any 



symmetric tensor A ab ym can be expressed as (see Section 5 of [8]) 



A 



ab 
sym 



Ag c 



(24) 



lB 



for some scalar A. L 3 ! The background Ricci tensor \jZ ah ^ ^ is of this form. Then, given that 



gbdK and TZ ad = -g bc K abcd , the form of Eq. (|23J) follows 



J^abcdef 



(0) 



is associated with a metric tensor. 



And so we have found that any index on 
Then, any index on 

V a hbc^ dhef 

must be contracted with one of the other five indices. 

Let us now choose the transverse and traceless gauge for the gravitons, B 



(25) 



Va^, h\ = 



(26) 



3 This statement is valid on the horizon's bifurcation surface. However, as Wald's integral expression can 

be evaluated over an arbitrarily chosen cross-sectional slice 6 , one can always calculate on the bifurcating 

slice without loss of generality. 

4 This choice is for convenience only, as the Wald entropy is gauge invariant. 
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So that contractions such as 



and 



A 



abed 



1 (0) _ _ 

Veh e a V b h cd 



A 



abed 



(0). 



V e h ab V c h e d 



(27) 



(2f 



vanish. 

To show that the term ( 1281) vanishes, we first argue that it must be of the form 



A 



bd 



(0) 



V e h f b V f h e d , 



otherwise a trace h a a appears. Integrating by parts with V e , we have 



A 



id 



(0) 



V e h f b V f h e 



A 



i,d 



(0) 



h f b V f h e d - 



A 



bd 



(0) 



h f h V e V f h e d 



(29) 



(30) 



All derivatives can be eliminated from the second term: 



V e V f h e d 



V f V e h e d + [V e ,V/]/l e 



PV«] (0) a 8 * + [V] IUJ ^ 



al (0) 



(31) 



which follows from Eq. ( 126|) and the standard identities relating commutators of derivatives 
to the Riemann tensor. 

Integrating by parts V/ in the first term in Eq. (!30]) : we find 



- V P 



.A 



6</ 



(0)\ . _ 

h f b V f h e d = V f V e 



A 



id 



1(0) 



h f b h e d + V e 



A 



bd 



V/ZiO /. e d , (32) 



so 



A' 



6(/ 



1(0) 



= [JB /e 6 f 0) h f b h* d ; 



(33) 



A 1 



bd 



where, besides transversality, we have redefined the background tensor [i3^ eM ] °^ = V' 
to emphasize that the above is really a mass term. Hence the term (1281) is a mass term. 

One can now see that the two V's in a kinetic term must be contracted, leading to the 
claimed form in Eq. (J2~T|) . But what about h 2 terms with greater numbers of derivatives? 
Terms having more than two V's can also make kinetic contributions, as follows. 



(0) 
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Derivatives have to be contracted in pairs since they cannot be contracted with a graviton 
index because of the gauge condition. Consequently, given a generic term with exactly 2n 
derivatives 



& 



a 1 a2...a2n-ia.2n 



(0) _ 

V Ql Va 2 • • • Va 2n _i Va 2 „ 



A 



abed 



(0) 



h„hh 



ab'^cd 



(34) 



one can use integration by parts, commutator relations like 



[V a , V 6 ] h cd = [K 



1(0) 



1(0) 



(35) 



symmetry properties of the background such as 



abcdl (°) 



[n abcd ] 



_ j^&ocdj (°) 
_ j^cotaj (°) 



_ j^oftdcj (0) 
^fecadj (°) 



defeal (°) 



[n dcba ] 



(36) 
(37) 



and the gauge conditions (1261) . to manipulate the expression into a power series in □ = V e V e 
acting on the gravitons. (Similar to our manipulation of Eq. (128j) into a series terminating at 
□°.) That is, 

n 

Y,[-A ahcd ]fh ab ^h cd . (38) 

3=0 

For instance, a term initially containing four derivatives and two gravitons can, after 
sufficient manipulations, be reduced to a combination of three forms: 



A 



abed 



1(0) 



h ab D 2 h cd , [A ahcd ]f h ab Uh cd , [A ahc X' KbKd . 



abcdl (°) 



(39) 



This outcome follows from our earlier discussion, where it was shown that the background 
geometry can only act on gravitons so as to contract indices or yield pairs of contracted 
derivatives. The essential point is that higher-derivative terms can contribute to the kinetic 
terms and therefore to the Wald entropy, provided that all but two of the derivatives act on 
the background tensors. 
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3.2 Expanding the Lagrangian 



Let us now begin the formal calculation by writing down the second-order expansion of ^J—gC: 



St® 



d{^C) 5g 



ub 



dg 



ab 



(2) 



+ 



+ 



dC 



on 



-6K, 



abed 



abed 



(2) 



+ 



dC 



1(2) 



d[V ai TZ a bcd. 



■oV a^abed 



dC 



where 



5C 



<9[V 

(ai V a,2)R'abcd\ 

_ S[y/=gC] 



'(ai V ' a 2 )Ttabcd 



(2) 



+ . . . , (40) 



(41) 



and the ellipsis has now be used to denote variations with respect to ever-increasing numbers 
of symmetrized derivatives. 

To proceed, we follow an iterative procedure that was first laid out in the third section of 
jl]. The basic idea is that a term like 



dC 



_<9[V( ai . . . V aj )T^abcd 

can always be re-expressed as 

dC 



(ai ■ ■ ■ » aj)T^abcd 



(2) 



1(2) 



V(ai 5V a2 . . . V aj )TZ a bcd 



(42) 



(43) 



(ai • • • V a j ) T^abcdj 

plus terms that are proportional to V ai Sg. Then, integrating everything by parts, one has 
(up to surface terms) 



dC 



l(2) 



5V a2 . . . V aj )T^abcd 



(44) 



(ai \d[V (ai . . .V aj) K abcd ] 
plus terms that are proportional to Sg. One can repeat this process j — 1 more times until 

(2) 



obtaining 



V(a. ■ ■ ■ V ai) 



dC 



8K, 



abed 



,^[V(ai ■ • • V ' aj )1Zabcd 

along with a collection of terms that are proportional to 5g (as well as surface terms). 
Consequently, we can reorganize the expansion (j40|) as follows (up to surface terms): 

Sim = [ W <* 6 g ab ]W + [ X ^8n abcd ] {2) , 



(45) 



(46) 
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where W a6 is some tensorial function of the geometry (with its precise form being irrelevant 

■ft 



to what follows □) and 



Mated \d[v ai n 



"abcd\ 



+ v '" v -' G[v ( J!,^ m1 ) + - < 47) 

is a tensor with the same symmetry properties of the Riemann tensor. Notice that X ahcd is 
the generalized version of the same-named tensor in Eq. (114]) . 

3.3 Isolating the kinetic terms 

Now, since 5g a b = g a b — 9 a °b > it follows that the complete non-expanded form of the W term 
in Eq. (|46p contains a factor g a b. Such a term can be dismissed, as integration by parts can 
be used to kill off any would-be kinetic contribution. [] 

To better understand this argument, suppose that we had the generic form 

[V (ai • • • V^h*] A***™*^ + h cd ) . (48) 

Here, all expressions are regarded as full unexpanded expressions except where indicated. To 
create a kinetic term, it is then necessary to move derivatives (via integration by parts) until 
obtaining a combination of the form 

V aj h ab [V^ • • • v a2 A abcd ^~ a i] V fll (^ } + h cd ) . (49) 

But such a term contains the vanishing factor V ai g c d = , and so the would-be kinetic term 
never has a chance to materialize. 



5 But note that, for C = C[g ab ,TZ abcd ] theories, W ab = dC/dg ab + \g ah C . 

6 Let us emphasize that this is different than starting with, say, STZ^ ~ VhVh + • • • and then integrating 
by parts to come up with /iVV/i . In this example, there may be an undifferentiated h, but it did not 
originate from a metric. 
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Having dismissed the metric variation, let us next focus on the 51Z abcd contribution. 
Schematically, this goes as 51Z = V5T + 5T5T , with ST = Vh . More precisely, 

<*ra (1) = \\v b h a c +v c h a b -v a h bc ] , (so) 



from which one obtains 

SK^Jh] = X - \V c V b h ad + W d V a h bc - V d V b h ac - V c W a h bd ] (51) 

and 

6n{ aL = \ \Vch ea V d h e b + V c h ea V h h° d - V e h ca Th db - V a h ce W d h e b 

-V a h ce V b h e d - V c h ea V e h db + V e h ca V d h e b + V e h ca V b h e d + V a h ce V e h db ] 

- {a — >d} . (52) 
From these expansions, it can be deduced that 

sl m = ^ (o) 6n m d + [xabcd] (i) 5n w d (53) 
= 1 [x***\ (0) ( V c h ea V d h\ + ■■■) + \ \x abcd ~\ (1) (W c W b h ad + •••), 



where the subscript k indicates that we only intend to retain the kinetic contributions and 
the ellipses denote the various permutations of the displayed indices. The first term on the 
right-hand side can clearly be kinetic and is what would normally be attributed to the Wald 
entropy. The question then is what becomes of the second term? 

To make sense of the second term, it is necessary to expand the tensor X ahcd to first order. 
Following the same iterative procedure as before, we have 



abed 
ypqrs 



^abcdj(l) = [ Z abcd;eq (0) ^(1) + 

where Z abcd ' e f is a tensor akin to W ab and 

abed dX abcd ( dX abcd 

ypqrs = _ _ Vfli 



(Q^abcd 
W^M ] + - (55) 
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is a tensor that is "two-fold" Riemann symmetric. 

Since Sg e f = g e f — gf) , we can, as previously discussed, neglect the first term on the 
right-hand side of Eq. fl54|) . Also recalling Eq. (l5"Tj) for 5TZ^\ we have 



abcdl W 



[x abcd ] 



abed 
ypqrs 



(0) 



(56) 



The insertion of Eq. f )56|) into Eq. (154]) then yields 



5£ 



(2) 



1^(0) (v Aa v^ + ...) 



+ 



abed 
ypqrs 



(0) 



(VcVfc/iad v r v, v + ■ ■ • ) • ( 57 ) 

As already mentioned, a four-V term might still make a kinetic contribution, so that the 
second term can not be dismissed. Nevertheless, we will now proceed to demonstrate that 
such a term can not contribute to the Wald entropy, simply because this is actually a mass 
term. 

In the Wald-entropy prescription [3j HJ E] the relevant metric perturbations are from the 
restricted class h a b\ n . ^ a b } = { rt } and 

h a b = V a Xb + VfeXa , (58) 

where Xa is a vector field that (at least) limits to the Killing vector on the horizon. Actually, 
as mentioned in Section 1, the Killing vector has been normalized such that 

tab = V a Xfe , (59) 

is the horizon binormal, for which 

e ab ^ iff {a^b} = {r,t} ; (60) 

and so only the off-diagonal elements h r t are relevant. 

Let us now consider a term such as V c VfeV a Xrf • O ne can use the Killing identity (Eq. (|E|)) 
to rewrite this as 



V c V b V a Xd = V c [n d abeX e ] 
= Ttdabe^ ' cX C ) 



(61) 
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where the second line is a consequence of \ a = on t ne horizon. This outcome leads us to, 
for instance, 

V c v b h ad = nf abe V cX e + nf dbe w cX e 

= rc£L Wcx e + v e Xc ] 

= K$LK, (62) 

where Eq. ( 1581) has been used twice and the second line comes about from the anti-symmetries 
of both the Riemann tensor and the horizon binormal vector; cf, Eq. (fTU]) . 

It now follows that the term proportional to y in Eq. (1571) can equivalently be written as 

(0) 



abed 
ypqrs 



(Kf a \ e KK^ qw K + •••); (63) 



so that any such term is a mass term as promised. Y_ 

Hence, the kinetic contribution and the Wald entropy is determined strictly by the first 
part of Eq. (|57|) . Following the result from Eq. (23) of [5] (which follows from Eq. (|52|) and 
the symmetry properties of [A' abc<i ] ^ ) 

[X ahcd ] (0) 5TZ^ cd = ~ [X abcd ] (0) (y e h bc V e h ad + 2Th ac V b h de ) (64) 

and recognizing that the second term can be gauged away, we arrive at 

5C { ? = \ \x ahcd f ] V e h bc V e h ad . (65) 

We can be even more precise by remembering that the relevant perturbations are sourced 
strictly by the off-diagonal elements in the {r, t} sector of gravitons. Hence, 

5tf = i[^ r4 ] (0) fyh.&.h*. (66) 



7 If the purpose is to calculate the Wald entropy, then one cannot "cheat" by having the Killing identity 
preceded by an integration by parts. That is, one is not permitted to turn SIZ^ ~ VhVh into hWh ~ hTZh 
and then argue that this is a mass term. 
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3.4 Interpretations 

Our claim is that all information about the Wald entropy is (up to normalization u) encoded 



in the horizon value of the single tensorial component [X 



sw — C [X abcd ~\ *■ ^ e a be c d 



rt l(0) 
rti ' 



(67) 



where sy/ is the "Wald entropy density" (or entropy per unit of horizon cross-sectional area), 
C is a universal normalization constant and we have made use of Eq. flBDl) to express the result 
in terms of the horizon binormal vectors. 

To compare, let us recall that the Wald entropy goes as 



swdA 



(6£ 



H 



where dA denotes an area element for a cross-section of the horizon "H. The density sw goes 
(in our notation) as 



s w = -2tt [X \ e ah e. 



cd ■ 



(69) 



And so the Wald formula agrees with our expression (IBTj) . with the normalization now fixed 
at C — —2n . 

Let us further clarify the relation between the Wald entropy and the coefficients of the 
kinetic graviton terms. We recall that, for the derivations of the Wald entropy [31 HJ [6], the 
idea was to start with the linearized field equation 



-g 



dg 



<ib 



h 



ab 







(70) 



and reduce this to a boundary term over a cross-section of the horizon. Following this path, 
one ends up with various terms of the generic form 



[A a ^ ab ] {0) V ai ...V aj h ab . 



(71) 



3 The correct normalization can always be uniquely fixed by the Einstein Lagrangian. In this case, X ah cd 



5 [9c9 b d ~ 9d9c] . and so [X 



rt l(0) 
rti 
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We now recall from the previous subsection that the relevant gravitons can, when pulled back 
to the horizon, be exchanged for e^'s or V a Xt 's (see Eqs. (I58f59p ). Meaning that the Killing 
identity (jSJ) can be used to reduce the number of derivatives when j > and, since the 
Killing vector vanishes on the horizon, any of these terms can be cast into the form 



Let us next consider the process of going from the original volume integral to a surface 
integral over the horizon and, subsequently, to a closed integral over a horizon cross-section. 
One is then required to apply Gauss' theorem twice. Starting with the integrated form of 
Eq. (1721) and following this route backwards, we have 



with A being the affine parameter for the horizon and M. the exterior spacetime. Although the 
first equality is trivial, the second equality is quite complicated, as sensitive limiting proce- 
dures are required to translate coordinates and geometric quantities from the null horizon to a 
timelike "stretched horizon" . The final form also assumes that the graviton and background 
do not depend on the "non-radial" spatial coordinates Xi, x^ ■ ■ -Xd-2- Yet, the underlying 
message is clear: Any contribution to the Wald entropy must necessarily come about from 
terms in the linearized field equation (or, equivalently, the quadratic action) carrying two 
derivatives. Our analysis explicitly establishes this connection and also identifies the origin 
of the kinetic terms. 

As observed elsewhere IS] , the Wald entropy can be generalized to other types of grav- 
itational couplings by a different choice of polarization for the gravitons. For instance, the 
shear viscosity of a black brane is determined by the kinetic coefficient of the h xy gravitons, 
where x and y are transverse directions on the brane that are mutually orthogonal as well as 
orthogonal to the direction of propagation. On this basis, it had been conjectured [9] that 



(72) 




(73) 
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the shear viscosity 77 could be determined in analogous fashion to the entropy; that is (cf, 



where e a b is a suitably defined binormal vector. 

However, our previous use of the Killing vector is unique to the h rt gravitons on the horizon 
and, hence, unique to the Wald entropy. For this reason, the calculation of any other type of 
coupling (such as 77) would generally involve the y term in Eq. (|57p . This term only becomes 
relevant for six (or higher) derivative theories. This is because, for a two (four) derivative 
theory, y = ( y ~ gg ); and the would-be kinetic terms are either identically zero or 
effectively zero through integration by parts. 

3.5 Summary 

In this section we have verified our assertion that, for a theory of gravity with any number of 
derivatives, the kinetic terms for the h rt gravitons completely account for the Wald entropy. 
By expanding out the Lagrangian to second order in gravitational perturbations, we have 
established that, due to the Killing identity, the only contributing terms are those for which 
a single component of the Riemann tensor is responsible for both of the gravitons. This 
property is essential to the applicability of the Wald formula to higher- derivative gravitational 
theories. Indeed, analogue formulas for other types of graviton coupling would, as discussed 
above, break down at the six- derivative order. 

4 Including matter 

Our considerations have so far been limited to theories of gravity without matter. It is then 
natural to ask if the inclusion of matter fields could alter any of the results of the preceding 
section. We now address this question and demonstrate that, even for a general theory of 
gravity coupled to matter, all our previous conclusions remain valid. 



Eq. (J67D) 




(74) 



a+b 
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4.1 Preliminaries 



Let us now add matter fields, denoted collectively by ip. Since anti-symmetric combinations 
of derivatives can always be replaced by 7^'s, the Lagrangian is of the form 

£ = £ [9ab, Kabcd, V ai 7^ af)cd , V( Ql V a2 )TZ abcd , . . . ; ij), V ol -0, V( ai V^)^, • • ■ ] ■ (75) 

In principle, one should then add the following set of terms to the expansion of in 
Eq. (@0D: 





(2) r 


dip 


+ 



dC 



Yd[V ai iP] 



SV ai iP 



(2) 



+ 



dC 



r5V (ai v a2) v> 



(2) 



+ 



(76) 



[d[V (ai V a2) iP] 

However, for scalar (vector, tensor) matter fields, at least the first four (three, two) terms 
in the series will not contribute. To understand this point, let us consider the case of three 
derivatives or less acting on a scalar <p. (Two derivatives acting on a vector and one, on a 
tensor would be equivalent situations.) As we have seen, a kinetic graviton term can only 
arise out of a variation of the Riemann tensor, 
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V5T + 5T5T where 5T ~ Vh 



(77) 



or out of variations of Christoffel symbols. But recall that only the 5T5T part of the Riemann 
variation contributes to the Wald entropy. Hence, a minimal requirement is having two 
covariant derivatives that act as Christoffel symbols. 

Now, for a scalar field, three-derivative terms can be dismissed because, as already stressed 
in Subsection 2.1, derivatives can only emerge out of the background in pairs. The presence 
of a scalar field (unlike a vector and other odd-spin fields) cannot viably alter this outcome. 
This leaves the remaining possibility of 



-4 ab V (a V, 



(78 



At a first glance, this seems to satisfy the minimal requirement. But, as the right-most 
derivative is required to act directly on 0, the above expression reduces to 



A ab V {a d b) <p 
19 



(79) 



and only one Christoffel symbol is available. It follows that such a term can only induce the 
variations 5g ab and 8<p\ thus disqualifying it as a kinetic contributor. 

However, the story could change if the matter sector had a sufficiently large number of 
symmetrized derivatives. The minimum requirements being at least four derivatives with a 
scalar field, three with a vector or two with a tensor. 

4.2 Four derivative terms 

Let us first consider the simplest example of possible contributions: 

= 0rt cd V (a V 6 V c V d) , (80) 

with 0, again, a scalar matter field. [] The variation of with respect to V( 0l V a2 V a3 V a4 )0 
will lead to a "candidate" kinetic term of the form 

[8C$ = [<Pg ah g cd f ] [5V (a V fe V c V d) 0] (2) . (81) 

A kinetic term might appear if any two of the derivatives act as a Christoffel symbol 
and the remaining two act directly on the scalar since such a combination leads to the the 
schematic form [5 (rTVV0)] (2) , and so T[h]T[h]W<j) ~ W(fVhVh . Nonetheless, explicit 
calculations have indicated that the net kinetic term from this Lagrangian vanishes. This 
finding can be explained by the following observations: Even though appears to have 
4! = 24 distinct terms, these can be joined into three types 

£ = 80D 2 + 80V a V fe V a V 6 + 80V a DV a . (82) 

One can then use the commutation relations (e.g., [V a , Vb]V c = Ti abc d V d ) to iteratively 
convert the second and third type into the first. For instance, one of the 24 terms goes as 

<pg ab g cd V a V c V b V d <p = <Pg ab g cd V a V b V c V d <p + <pg ab g cd V a [V c ,V b ]V d <p 

= 0D 2 + <Pg ab g cd V a n cbd e V^; (83) 
9 The second scalar field in front is to prevent this Lagrangian from being a total derivative. 
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and similarly for the other 15 terms that are not initially of the </>D 2 (/> type. 
The end result of the just- described process is 

£ = 240D 2 + <Pg {ab g cd \ a {Tl bcd e V e 0) , (84) 

with symmetrized indices on the metrics. Although there (again) appears to be 4! = 24 
different terms in the Riemann part of the above expression, this really only contains 16 such 
terms. The reason being that 8 of the 24 terms are of the form 

4>g ad g hc V a K hcd e V e , (85) 

which is already identically vanishing through the contraction of the first two Riemann indices. 

The rest of the Riemann part vanishes identically as well. To understand this, let us 
suppose that the index on the derivative V a is fixed while the other 3 Riemann indices bed 
remain symmetrized. This leads to the following 6 terms (where we display only the Riemann 
tensor for brevity): 

Tt(bcd) — T^-bcd e + T^-dbc e + Tt-cdb e + Tt-bdc e + Ttcbd e + Ttdcb e ■ (86) 

But, the first three terms sum to zero and, likewise, the latter three terms, due to the Jacobi 
identity; cf, Eq.( l37|) . Then, since the sum total is (as we vary the index on V) four such 
vanishing sets, the Riemann part of Eq. (184")) is zero. 

What is left to establish is that [(Jn 2 ^]^ similarly makes no kinetic contribution. This 
can, indeed, be verified with an explicit calculation but also follows from a simple argument. 
To show this, let us first recall that 

1 

^9 

After we disregard all kinetic contributions that can be gauged away, which include either 
derivatives acting on determinants or derivatives acting transversely, one finds that the sole 
potential contributor is 

□ 2 -> g ab [d a d b g cd ] [d c d d( j>] . (88) 
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□ 2 = ~^d a {g ab V=gd b 

v^g I 



--d c {g cd ^gd d(l 



7) 



However, even this term can still be gauged away through integration by parts, and so van- 
ishes. 

4.3 2n derivative terms 

The same trend extends to any number of symmetrized derivatives. Again working with a 
scalar matter field for simplicity, we can justify this claim via the following arguments: 
Let us begin by considering the Lagrangian 

C 2 M = <j>g a ^ ■ ■ ■ g a ^'^V {ai V a2 . . . V a2n _ x V a2n) , (89) 

which contains 2n symmetrized derivatives. Similarly to the four-derivative case, £2 n {4>) 
should decompose into the form 

n-l 

0LT0 + [n [k] V [2n - 2k - 1] ] a V a ct) , (90) 

k=l 

with the square bracket meant to represent a collection of k (4-index) Riemann tensors fol- 
lowed by 2n — 2k — 1 V's contracted in all possible ways. Such an arrangement is always 
possible given a sufficient number of a commutations of derivatives, 

Now, the key point is that any of the generated Riemann tensors arises due to a commu- 
tation of symmetrized V's, and so is of the basic form 

[v (a4 , v a ,] v 0fc) . . . = n {aiajak) e v e . . . , (9i) 

whereby three of the Riemann indices are explicitly symmetrized while the fourth index is 
summed over independently. So that, just like for the 4-derivative example (Eq. (I86I) ) any 
Riemann tensor produces two sets of 3 terms with each vanishing due to the Jacobi identity. 

Let us now focus on the 0D n term. Again disregarding kinetic contributions that can 
be gauged away, we are left with potential contributors of only the two basic forms, 

. . . g ab g cd d a h^d c h^ . . . dbd&dfdidrf and ... g ab d a h cd d b W . . . d c d d d,d 3 <P , (92) 

10 To ensure the displayed ordering, one should commute derivatives from left to right; i.e., opposite to the 
direction of Eq. (|83|) . 
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with the dots indicating other irrelevant structure. The crucial point here is that the deriva- 
tives act only symmetrically on the gravitons, whereas a Riemann tensor is constructed out of 
anti-symmetric combinations of derivatives. Meaning that these terms are simply not capable 
of hiding a Riemann tensor. 

4.4 Summary 

Given a generalized theory of gravity, we have now shown that adding matter fields with any 
number of symmetrized covariant derivatives acting on them does not lead to any kinetic 
contributions beyond those already encountered in Section 3. The explicit calculations used 
scalar fields but can, as discussed, be generalized in a straightforward manner to any type of 
tensor field. 

5 The generalized field equation 

The calculation of the quadratic action for a generalized theory is often complicated. But one 
can rather use the linearized field equation as an equivalent but simpler means for extracting 
the Wald entropy. The compatibility of the two approaches follows from their equivalency 
up to total derivatives. An appropriate choice of surface term can always be used to cancel a 
total derivative and then, as explained in [6], any such boundary term does not contribute to 
the Wald formula. 

The entropy should be extracted from terms of the schematic form VV/i in the field 
equation, subject to the various subtleties discussed. The field equation has already been 
presented (without a full derivation) in [1]. However, for practical calculations, we have 
identified a form of the equation differing by a sign from that of jl]. 



23 



5.1 The linearized field equation 



Let us presume, for simplicity in the presentation, that the matter fields ip are coupled only to 
the metric and any such terms have been collected separately in \C^. Then the gravitational 
field equation for the density 



C 



total 



-fj 



£ + 2 A* 



is 



where 



<9£ 1 

— 5g pq + -g pq C5g pq + X abcd 5K abcd [5g pq ] + Q pq 5g 



-T pq 5g. 



VI i 



rppq ^ 



2 d[y/=gC^] 



(93) 



(94) 



(95) 



V9 dg pq 

X abcd is defined in Eq. fj47j) and the tensor Q pq accounts for the extraneous terms in the 
iterative procedure of Section 2. (But, if C contains only anti-symmetrized combinations of 
derivatives, there would be no Q contribution.) Note that, had we varied with respect to the 
contravariant form of the metric Sg pq , then the stress tensor needs to be defined as 

2 dW=gC^[ 



T — 

J-pq — 



V9 dgPi 
Now linearizing and rearranging, we have 

dC 



(96) 



abcdl (°) 



[x abcd ] 



J'VQ _|_ 



1 



-g pq C + Q 



(0) 



b VQ i 



(97) 



dgpq 2- 

with the right-hand side being irrelevant to the Wald entropy. 

We will thus focus on the left-hand side and recall the expansion for STZ^ cd in Eq. ( |5~T|) . 
This expression and the symmetry properties of the background Riemann tensor (which are 
shared by [A' a6cd ] ^° ) indicate that the left-hand side of Eq. (19~T|) reduces to the sum of four 
equivalent terms. Denoting this sum as Q, we obtain 



[X abcd } {0) 5rt 1] 



2[X apqb f ) V a V b h pg 

g pq h P q . 



abed l lb Vq 



(9* 
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For the rest of this section, it is implied that the zeroth-order geometry applies to all tensors 
besides the h's. 

To proceed, one considers separately the symmetric |{V a ,Vb} and the anti-symmetric 
|[V a ,Vfe] combinations of the derivatives. Whereas the former is trivially handled with a 
double integration by parts, the latter is more complicated. Nonetheless, it is possible to con- 
vert {V a , Vb} into Riemann tensors. Then, through repeated application of the symmetries 
of both the X and 1Z tensors, one obtains the expression in the second line of Eq. (|98|) . What 
is left is to strip off the graviton and then linearize Q pq . 

Skipping over subtleties, we find that the complete field equation now goes as (cf, Eq. (1971) ) 

2V b V a X apqb - X abcp TZ abc q + -g pq £ = -T pq , (99) 

where we have assumed that £ = £ [TZ a bcd\ for simplicity. This is in contrast to the relative 
positive sign between the first terms in Eq. f )99|) in [I]. 

5.2 Simple examples 

To clarify this process, let us recall the Einstein Lagrangian, C-Ein = 1Z , regarded as 
independent of the metric. Plugging this into Eq. ( HTj) . we find that 

= \ [9 aC 9 bd ~ 9 ad 9 bC ] , (100) 
and the tensor Q pq of Eq. ( ]9~8]) then reduces to 

Q P E \n = ~K Pq ■ (101) 

This is indeed the correct form, as substituting back into the field equation (|99|) . we obtain 

K pq - -g pq Tl = T pq . (102) 

Let us further confirm the consistency of Eq. ( 198]) by starting with an £ = F(TZ) theory 
of gravity ( 1Z = g ac g bd, R-abcd ) and following the standard procedure, 

d n F{K) [n pq 5g pq + g pq SK m ] - ±F(n)g pg 5g pq = T pg 5g pq . (103) 
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To handle the Ricci variation, one can use a contracted form of Eq. ( 15TT) and then integrate 
by parts (twice) to free up the graviton. In this way, one ultimately finds that (e.g., [TO] ) 



or, equivalently, 



[K pq + g pq D - V p V g ] d n F{K) - \g pq T(Jl) = +T pq (104) 



[V^V 9 - g pq D - K pq ) d n F{TZ) + ^g vq T(Jl) = —T pq . (105) 



One can verify that this is in perfect agreement with Eq. fl99|) . given that 



p^apqb 



dF{TZ) 

l -[g aq g pb -g ab g pq ]d n F{K). (106) 



6 Example calculations 

6.1 General considerations 

We wish to demonstrate by specific examples that the coefficients of the h rt kinetic terms can 
indeed be used to directly extract the Wald entropy. We will use two gravitational theories, 

C Q = K + aU ahcd n abcd , (107) 

Cp = n + (3V k R abcd V k R abcd , (108) 

where a and are constants. From a physically motivated perspective, these parameters 
should be regarded as small: a Cr^ and <C r\ , where is the horizon radius of the 
corresponding black hole solution. 

The Einstein term TZ allows us to normalize our results and, more importantly, given that 
the corrections to Einstein gravity are small, it allows us to use the Einstein background 
solution in the calculations. Any contribution from a perturbative correction is already first 
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order in a or /3, so that the Einstein background solution suffices. The kinetic coefficients in 
Einstein's theory are numerical constants and so are insensitive to the form of the solution. |^ 
In the following, we incorporate the notation g ab = , TZ abcd = H abcd and r «&,c = , 
where T ab:C = g cd T d ab . 



6.2 4-derivative gravity 

Let us first consider the theory Eq. (1107!) . Keeping in mind that the Lagrangian is integrated 
to get an action, I = J y/^gC a d D x, we can integrate by parts. 

To begin, we write the linear in a and quadratic in h part of the Lagrangian density 

5[^gC a f ] = aS[^gg aa g bl3 g^g dS n a(3j5 n abcd \ {2) (109) 

and look for any term that makes a kinetic contribution. As a first step, let us consider the 
"(1,1)" terms or 5 [y/^gg aa g b ^ g^g^TZa^s] oTZ^^. Our previous analysis suggests that it 
does not contribute. We will demonstrate this explicitly by working with the unexpanded 
version of the factor \[--gg aa g h ^ g n g^'R-ap^- So we rewrite Eq. fl 109[) as 

5 [^/g£ a } {2) = = aT^VVV'X/^ \V b V c h ad - V b V d h ac - V a V c h bd + V a V d h bc ] , 

(110) 

where Eq. floT]) has been used to expand out STZ^ cd . 

Now, to order h in the square brackets, any of the V's can be replaced with a "full" V. 
Let us do so with the left-most V in each of the four terms and then integrate by parts: 

- a^gg aa g hfi g n g d& [V b K a p y s^ ' c h ad - V^^V ' d h ac - W a H a p 1& W c h bd + V ' aTZ a p 1& V ' d h bc ] . 

(Ill) 

As in Section 2, when two V's act on a graviton, the Killing relation can be used to 

convert them to a Riemann tensor. So, to get a kinetic term out of Eq. (II lip , we need to 

n Here, it is assumed that the intention is to calculate the Wald entropy in units of horizon area. The 
Einstein contribution to the horizon area does get corrected. 
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expand TZ a /Sj8 to first order in h and to zeroth order in derivatives. Then, since 7Z a bcd = 
db^ ac ,d — d a Tbc,d + g mn ^ ac,uX nb,d — g mn ^bc,mTna,d and r a ;, )C = | {d a gb c + dbg ac — d c g a b) , the 
only terms in TZ a b c d that can possibly contribute are g mn T ac ^ m T n b td and g mn Tb c ,mFna,d ■ But 
Vbg mn = V a g mn = , so there can be no kinetic contribution at all. 
Hence, we have 

6 [V££a] (2) = 2a [V^n abcd ] (0) 5U^ cd + ■ ■ • , (112) 

or, using Eq. f lMl) . 

5[^gC a } {2) = 2a^n abcd ^ m h bc V m h ad + 2V m h ac V b h dm ] , (113) 

where the second term in the square brackets can be gauged away. 

We can identify sw = 2aC1Z abcd e a &e C( j , with the normalization C = —2ir given by the 
Einstein term. This is in agreement with that obtained via a direct application of Wald's 
formula. 

6.3 6-derivative gravity 

The Lagrangian is now given by Eq. (11080 . Here, we begin the quadratic density 

5[^C a f ] = (35[^g aa g b Pg^g d5 V k K a ^ k K a b cd ] {2) . (114) 

The (1,1) terms, again, cannot contribute. Because of the derivatives in front of the Riemann 
tensors, either of these must be expanded to first order in h and zeroth order in derivatives. 
Hence, by the very same reasoning as provided above, a kinetic contribution can not be 
obtained. 

This leaves us, after integration by parts, with 

5 [^] {2) = -2/3 [v^V fc V fc ^] (0) SR^d + ■ ■ ■ » ( 115 ) 

or, like before, 

6 \^C a f ] = -2/3v/^V fc V fc 7T 6cd \W m h bc V m h ad + 2V m h ac V b h dm ] + ■ ■ • , (116) 
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and so the identification s\y = —2/30071 tabbed follows. This is, once again, in agreement 
with the Wald formula when the normalization is C = —2n . 

Finally, more derivatives beyond six would neither conceptually nor technically complicate 
the calculation. 

7 Conclusion 

To summarize, we have investigated Wald's Noether charge entropy [H 0], relying on its 
identification with a quarter of the horizon area in units of the effective gravitational coupling, 
as first established in [5]. The Wald entropy can, as now verified, be determined from the 
kinetic coefficients for the h rt gravitons on the horizon. We have also clarified what terms in 
the quadratically expanded action (or linearized field equation) can contribute to the entropy 
and illustrated our general procedure with some of explicit examples. Additionally, we have 
reconsidered the gravitational field equation for a general theory of gravity and presented it 
in a form which differs from jl]. 
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